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Abstract 

We present an engineered version of the divide-and-conquer algorithm for finding the clos- 
est pair of points, within a given set of points in the XY-plane. In this version of the algorithm, 
only two pairwise comparisons are required in the combine step, for each point that lies in the 
25-wide vertical slab. The correctness of the algorithm is shown for all Minkowski distances 
with I. We also show empirically that, although the time complexity of the algorithm is 
still 0{n Ig n), the reduction in the total number of comparisons leads to a significant reduction 
in the total execution time, for inputs with size sufficiently large. 



1 Introduction 

The Closest-Pair problem is considered an "easy" Closest-Point problem, in the sense that there 
are a number of other geometric problems (e.g. nearest neighbors and minimal spanning trees) 
that find the closest pair as part of their solution fB entley[ [1980" p. 226]. This problem and its 
generalizations arise in areas such as statistics, pattern recognition and molecular biology. 

At present time, many algorithms are known for solving the Closest-Pair problem in any di- 



mension /: > 2, with optimal time complexity (see |Smid| [ |2000[ | for an overview of Closest-Point 
problem algorithms and generalizations). The Closest-Pair is also one of the first non-trivial com- 
putational problems that was solved efficiently using the divide-and-conquer strategy and it be- 
came since a classical, textbook example for this technique. 

In this paper we consider only algorithms for the Closest-Pair problem that can be imple- 
mented in the algebraic computation tree model. For this model any algorithm has time complex- 
ity Q.{n Ig n). With more powerful machine models, where randomization, ihQ floor function, and 



indirect addressing are available, faster algorithms can be designed |Smid[[2000l |. 
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Historical Background 



An algorithm with optimal time complexity 0{n Ig n) for solving the Closest- Pair problem in 
the planar case appeared for the first time in 1975, in a computational geometry classic paper by 
Shamos|fl975] . This algorithm was based on the Voronoi polygons. 

The first optimal algorithm for solving the Closest-Pair problem in any dimension /: ^ 2 is due 



to |Bentley and Shamos| | |1976| . Using a divide-and-conquer approach to initially solve the problem 
in the planej^ those authors were able to generalize the planar process to higher dimensions by 
exploring a sparsity condition induced over the set of points in the /:-plane. 

For the planar case, the original procedure and other versions of the divide-and-conquer algo- 
rithm usually compute at least seven pairwise comparisons for each point in the central slab, within 
the combine step (see |Kleinberg and Tardos| [ [2005p , |Bentley and Shamos|p976l l, and |Cormen et al. 
| [200T1 , for instance). 

In 1998, Zhou, Xiong, and Zhi]^ presented an improved version of the planar procedure, where 
at most four pairwise comparisons need to be considered in the combine step, for each point lying 
on the left side (alternatively, on the right side) of the central slab. In the same article, Zhou et al. 
introduced the "complexity of computing distances", which measures "the number of Euclidean 
distances to compute by a closest-pair algorithm" [ Jiang and Gillespiel|2007] . The core idea behind 
this definition is that, since the Euclidean distance is usually more expensive than other basic 
operations, it may be possible to achieve significant efficiency improvements by reducing this 
complexity measure. 

More recently, Ge, Wang, and Zhu used some sophisticated geometric arguments to show 
that it is always possible to discard one of the four pairwise comparisons in the combine step, 
thus reducing significantly the complexity of computing distances, and presented their enhanced 



version of the Closest-Pair algorithm, accordingly [Ge et al.[|2006 | . 

In 2007, Jiang and Gillespie presented another version of the Closest-Pair divide-and-conquer 
algorithm which reduced the complexity of computing distances by a logarithmic factor. However, 
after performing some algorithmic experimentation, the authors found that, albeit this reduction. 



the new algorithm was "the slowest among the four algorithms" ] Jiang and GillespieH2007] that 
were included in the comparative study. The experimental results also showed that the fastest 
among the four algorithms was in fact a procedure named Basic-2, where two pairwise compar- 
isons are required in the combine step, for each point that lies in the central slab and, therefore, has 
a relative high complexity of computing distances. The authors conclude that the simpler design in 
the combine step, and a consequent correct balance in trading expensive operations with cheaper 
ones are the main factors for explaining the success of the Basic-2 algorithm. 

In this paper we present a detailed version of the Basic-2 algorithm. We show that only two 
pairwise comparisons are required in the combine step, for each point that lies in the central slab, 
and that this number of comparisons is minimal. This result and the subsequent correctness of the 
Basic-2 algorithm is shown for all Minkowski distance^ with \. 



1 ^ 



According to Bentley |1980|, Shamos attributes the di scovery of this pr oce dure to H.R. Strong. 



Ge et al. 



|2006| and [Jiang and Gillespie| | |20Q7| for some 



^The article in question was published in Chinese. See 
expHcit references. 

^In fairness to all parts involved, we must say that all the main results presented in this paper, including the design 
of the Basic-2 algorithm, were obtained in a completely independent fashion from any previous work by other authors. 



It was only during the process of putting our ideas in to writing that we came across the articles by Zhou et al, Ge et al. 
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The rest of the paper is organized as follows. In Section [2[ we review the classic Closest- 
Pair algorithm as presented by |Bentley and Shamosmi976| |. In Section [3| we present our detailed 
version of the Basic-2 algorithm and give the correspondent proof of correctness. In Section |4} we 
present a comparative empirical study between the classic Closest-Pair and the Basic-2 algorithm, 
and discuss the experimental results obtained with distinct Minkowski distances. 



2 The divide-and-conquer algorithm in the plane 



The following algorithm for solving the planar version of the Closest-Pair problem was first pre- 



sented by [Bentley and Shamos] [ |1976 1 . 

Let P be a set of /I > 2 points in the XY-plane. The closest pair in P can be found in 0{n Ig n) 
time using the divide-and-conquer algorithm shown in Figure [TJ 



Closest-Pair(P) 

1 Presort points in P along the x-coordinate. 

2 Split the ordered set P into two equal-sized^ subsets by the vertical line /, 
defined by the equation x = Xmedian- 

3 Solve the problem recursively in the left and right subsets. This will give the 
left-side and right-side minimal distances and dji^ respectively. 

4 Find the minimal distance d^R among the pair of points in which one point 
lies on the left of the dividing vertical line and the other point lies to the right. 

5 The final answer is the minimum between d^^ dji^ and dui. 



"give or take (exactly) one point, when the number of points is odd. 



Figure 1 : Pseudocode for the divide-and-conquer Closest-Pair algorithm, first presented by [Bentley 
and Shamos|[[T976l . 



Since we are splitting a set of n points in two sets of n/2 points each, the recurrence relation 
describing the running time of the CLOSEST- Pair algorithm is T{n) = 2T(n/2) + f{n), where 
f{n) is the running time for finding the distance diR in step 4. 

At first sight it seems that something of the order of ^^/4 distance comparisons will be required 



to compute diR. However, Bentley and Shamos] [1976| | noted that the knowledge of both distances 
di and dR induces a sparsity condition over the set P. 

Let 5 = min{dL^ dR) and consider the vertical slab of width 25 centered at line /. If there is any 
pair in P closer than 5, both points of the pair must lie on opposite sides within the slab. Also, 
because the minimum separation distance of points on either side of / is 5, any square region of 
the s lab, with side 25, "can contain at most a constant number c of points" [ [Bentley and Shamos 
1976J , depending on the used metric]^ 



As a consequence of this sparsity condition, if the points in P are presorted by y-coordinate, 
the computation of diR can be done in linear time. Therefore, we obtain the recurrence relation 
T{n) =2T{n/2) + 0{n), giving an 0{n Ig n) asymptotically optimal algorithm. 



j2Q06[ , and [Jiang and Gillespie 
Tn the original article, Bent 



2007 



ey and 



which, obviou sly, take precedence and deserve due credit. 



Shamos |1976|, used the Minkowski distance doo and obtained the value c = 12. 
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Figure 2: "Hopscotching" A-points in ascending order. For each point visited on either side, 
the BASIC-2.S4 algorithm computes the distance for the two closer, but not lower, points on the 
opposite side. 

3 The Basic-2 algorithm 

In this Section we discuss a detailed version of the Basic-2 algorithm, which was first presented by 



Jiang and Gillespie [2007| 



The Basic-2 algorithm is an optimized version of the Bentley and Shamos procedure for the 
planar case discussed in section|2} In fact, the Basic-2 algorithm is the same as the CLOSEST- Pair 
algorithm (see Figure [T]), with the sole difference that the computation of the distance dui in step 
4 now requires only two pairwise comparisons per point to find the closest pair within the central 
slab. The pseudocode for computing the du^ distance in the Basic-2 algorithm is shown in Figure 

The time complexity of the BASIC-2.S4 algorithm is obviously 0{n), since it traverses once 
the arrays and F/^ in a "hopscotch" manner (see Figure [2]), and it takes constant time on each 
iteration. 

Since we are only interested in performing step 4 of the Closest-Pair algorithm, in the following 
we assume that x^nediam ^l, and dR are already computed. We also assume that the array Y contains 
a y-sorted partition of all points in P, i.e., the first and second halves of Y contain the points that 
are to the left and to the right of /, respectively, and both halves are sorted along the y-coordinate[^ 

We denote the vertical slab centered at line x = x^edian of width 25 by the symbol A and the 
central line by /. 

Before we show the correctness of the algorithm we first prove the following 



^The structure of the array Y may seem difficult to obtain, and therefore, be an extra source of complexity in the 
overall algorithm. However, this is not the case because the structure arises as a natural consequence from the need to 
maintain the y-presorting throughout the recursive calls in 0{n) time. 
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BASIC-2.S4(x^edzan, ^^L, Y) 


1 




— min{dL,dR) 


2 


Yl ^ set of y-sorted points in A lying to the left of / 


3 




set of y-sorted points in A lying to the right of / 


4 


left ^ first[YL] 


5 


right 


^ first[YR] 


6 


while there are points in Yl and Yr 


7 




do dist ^ distance{l eft, right) 


8 




if dist < dmin 


9 




then closestPair ^ {left, right) 


10 




dmin dist 


11 




if Vcoordileft] < Vcoordiright] 


12 




then if there is at least one more point in Yr 


13 




then dist ^ distance{left,next[YR]) 


14 




if dist < dmin 


15 




then closestPair ^ {left,next[YR]) 


16 




dmin dist 


17 




left ^ nextlYi] 


18 




else if there is at least one more point in Yl 


19 




then dist ^ distance{next[YL], right) 


20 




if dist < dmin 


21 




then closestPair ^ {next[YL],rigth) 


22 




dmin dist 


23 




right ^ next[YR] 


24 


return closestPair and dmin 



Figure 3: Pseudocode for step 4 of the Basic-2 algorithm. 



Lemma 1. Let rf^ : R x ^ R denote the Minkowski p-distance, I ^ p ^ oo. Let Pq = 
(•^05 Jo) ^ Yi{respectively^YR) be an arbitrary point lying in the central slab A, and let Fq ^ 
YR^respectively^Yi) be the array of points that lie opposite and above Pq, sorted along the y- 
coordinate. The closest point to Pq, in respect to dp, is either the first or the second element of 
Jo. 

Proof We first give proof for the Minkowski distance for = 1, defined by d\ (A,5) = \xa — ^^^I + 
\yA —JeV We assume, without loss of generality, that Pq = (O5O) ^ and, as a consequence, 
that Yo<ZYr. Let A — {a^b)^ ^ a ^ 5, Z? > be the first point in Fq- We note that, because Fq 
is sorted along the y-coordinate, it is sufficient to consider the case where £^ = 0, since all others 
cases with b>Q can be obtained by making an upper translation. This translation does not disrupts 
the relative positions of the elements in Yq and therefore, all arguments presented for = will 
remain valid. So, let A = {a^Q) and let P = (x,_y), ^ x ^ 5, _y > 0, be any other point in Fq- We 
must consider three cases (ilustrated in Figures |4} [5} and [6]). 
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Figure 4: Possible location for the first three points, A, P, and Q in Case 1. When A is centered it 
is possible for both A and Q to be the closest points to Pq, in Yq. However, this is a hmit-case. 



Case 1: a 



We find that 



^ 0-- 
2 



2 2 



6 5 

- ^ - 
2 2 

-l< - 
2' ^ 2 



On the other hand, we have 



Therefore, 



di{A,P)^5^ 



5 

2 



8 

2 



5 5 

di{Po,P)=x + y^x+-^ - = Ji(Po,A) 



which is to say that A is a closest point to Pq, in Yq. 

We note that, if we take J = f , the first three points in Yq may have coordinates 
A = (|,0); P = (5, |); 2 = (0, |), respectively. This is the limit-case depicted in Figure Q 
where A and Q are both the closest points to Pq, in Yq. 
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Figure 5: Possible location for the first three points, A, P, and Q in Case 2. When A is close to Pq 
all other points in Yq are "pushed" away by the sparsity of A. 

Case 2: ^ a < |. We consider two possibilities 



i) Let}; > |, then 



ii) Let J ^ |, then 



5 5 

di{Po,P) =x + y>x+-^ - >a = di{Po,A)' 



di{A,P) ^ 5 <^ \x-a\ 5 <^ 

I I / ^ / 2 2 



<^ X — <2 ^ — V X — ^2 ^ <^ 

2 2 

<;=^ x^a-\ — V x^a <0 = 

2 2 

^ V ' 

Contradiction 



x^ a + 



Therefore, 



5 5 

di{Po,P)=x + y^a+-+y^ ->a = di{Po,A). 



Considering i) and ii) we conclude that A is the closest point to Pq, in Yq. 
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Figure 6: Possible location for the first three points, A, P, and Q in Case 3. When A lies further 
from Po it is possible for another point to be closer to Pq- However, this point is necessarily the 
next lowest point, which is to say, this point is the second element in Fq- 



Case 3: | < a ^ 5. We must consider two possibilities 



i) Let a, then 

di{Po,P) =x + y ^a + y^a = di{Po,A). 
and, as in the previous cases, A is the closest point to Pq, in Fo- 



il) Let X < a, then 

di{A,P)^ 5 ^a-x + y^ 5 5-a+x 

This means that it is possible to have at least one point P = {x,y) G Yq such that 
di{Po^P) = X + y < a = d\{Po,A), and ^ x < a, 5 — a+x^y^a — x. Let 
2 = {x\,y\) G Fo, and assume that ji ^ }^ (i.e., Q precedes P in Fq). We know that 

di{A,Q) ^ 5 ^\xi-a\ +yi ^ 5 ^ 
<^ x\^a + 5—y\ V xi^a — 5+yi (1) 

From the first inequality of ([T]) we find that 

x\ ^ a + 5 —y\ ^a + 5— y>x + y + 5— y = 5+ x^5 

and this means that, in this case, for Q to precede P in Fq, it must lie outside the A slab, 
which is a contradiction. 

The second inequality of ([T]) holds only if we choose 5 — a ^yi ^y io guarantee that 
^ xi ^ yi — {5 — a). However, for this choice of the Q coordinates, and from the 



8 



sparsity of A we get 





di{P,Q)^5^\xi-x\+y- 


y\ 




<^ 


x\-x^ 5 +y\ -y 


V 


x\—x^ —5+y — y\ ^ 


<^ 


x\ ^ x + 5 + _yi — (5 — a+x) 


V 


^ {x + y)-5-yi ^ 


<^ 


x\^y\+a 


V 


xi < a — 5 —yi ^ 


<^ 


^\ > y\+a — 5 

^ V ' 


V 


xi < a — 5 ^0 



Contradiction 



The first inequality of (2) is a contradiction with our current working hypothesis (i.e., 
with the second inequahty of ([T])) and the second inequahty of (2) imphes that Q ^Yq. 
So, if there is one point P G Fq closer to Pq than point A, then no other point, but A, can 
precede P in Fq- Now, suppose that there is another point 2 G Fq that is also closer to Pq 
than point A. This means that, as with P, no other point can precede 2 in Fq- However, 
since both P and Q are in Fq, one must precede the other, which is a contradiction. 
Therefore, in this case, the only point possibly closer to Pq than point A is the second 
element of Fq. 

From all the previous three cases we may conclude that the closest point to Pq, in Fq, is either 
the first or the second element of Fq. This proves Lemma [T] for the Minkowski distance d\. 




Figure 7: Unit circles for various Minkowski p-distances. 



To obtain proof for all other Minkowski distances dp, p > 1, we take into account the fact that 
the convex neighborhoods generated by the Minkowski distances possess a very straightforward 
order relation, where larger values of p correspond to larger unit circles, as shown in FigurejT} This 
ordering means that the sparsity effect within the A slab will be similar, but somewhat stronger, for 
larger values of p. Therefore, the precedent analysis of d\ not only remains valid for > 1 but, 
in a sense, the corresponding geometric relations between elements of Fq are expected to be more 
"tight" for all other Minkowski distances. □ 

The preceding Lemma [T] establishes an upper bound of two for the number of pairwise com- 
parisons needed to be computed for each point in the A slab. Also, the analysis made in Case 3 of 
the corresponding proof shows that this bound is tight, i.e., we have established the following 

Corollary 2. For each point in the central slab A, the minimum number of pairwise comparisons 
required to compute diR, in the divide-and-conquer CLOSEST- PAIR algorithm, is two. 
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Correctness 

To prove the correctness of the BASIC-2.S4 algorithm we consider the following 

Loop Invariant. At the start of each iteration of the main loop, left and right are references to the 
points of Yl and Yr, respectively, with minimum y-coordinates, that still need to be compared with 
points on the opposite side. Also, dmin corresponds to the value of the minimum distance found 
among all pairs of points previously checked. 

We show that all three properties Initialization, Maintenance, and Termination hold for this 
loop invariant. 

Initialization At the start of the first iteration, the loop invariant holds since left and right are 
references to the first elements of Y^ and Yr, and both arrays are y-sorted in ascending order 
(by construction). Also, no pair with points on opposite sides of line / was checked yet, so 
the current minimum distance is the minimum between the left-side and right-side minimal 
distances dimin and dRmim respectively. This value is stored in dmin- 

Maintenance Assuming that the loop invariant holds for all previous iterations we now enter the 
next iteration. The first thing the loop does is to compute the distance between the points 
referenced by left and right and, if that is the case, it updates the value of dmin- Next, the loop 
determines which of the two points has smaller y-coordinate. Let us assume, without loss of 
generality, that in this iteration left is the lowest point. Since left G Fl, the algorithm checks 
to see if there is at least one more point in Yr (denoted by next\YR]), following the point 
right. If there is such point, the algorithm computes the distance between left and next\YR\, 
and updates the value of dmin accordingly. 

By our hypothesis we know that right and next\YR\ are the points on the right-hand side 
with the smallest y-coordinates that are still greater, at most equal, to leff^ y-coordinate. 
Therefore, and taking into account Lemma [TJ we conclude that the value of dmin corresponds 
to the minimum between the previous minimal distance and the minimum distance computed 
for all pairs of points that contain the left point. 

The iteration ends by incrementing the reference left to the next point in Y^, which means 
that the original left point will no longer be available for comparison. Note also that the right 
reference remains the same so that the corresponding point will still be compared with other 
points on future iterations (see Figure [8]). We may conclude that the new pair of references 
left and right, and the new value dmin still satisfy the loop invariant. 

Termination The loop ends when one of the references, left or right, reaches the end of the cor- 
responding array, Yi or Yr, respectively. Let us assume, without loss of generality, that the 
left reference reaches the end of the array Y^ and terminates the loop. This means that it was 
the left reference that was incremented at the last iteration and so, it was this reference that 
corresponded to the lowest point. Accordingly, the loop computed the distances between left 
and the two closer, but not lower, points in Yr and updated dmin- As a consequence, all re- 
maining pairs of points are composed by the point left and points that belong to the array Yr 
and lie in higher, more distant positions. Therefore, we have computed all distances between 
pairs of opposite points that may lie at a distance smaller than the current minimal distance, 
and so the value dmin corresponds to the minimal distance between all pairs of points in A. 
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Figure 8: Some iterations of the main loop in BASIC-2.S4. (a) Compute distance between left and 
right, the first elements of Yl and Yr, respectively, (b)-(c) Point left is lower than point right so, 
compute distance between left and the next element in Yr. Update left to the next element in Yl. 
Compute distance between left and right, (d)-(e) Point right is lower than point left so, compute 
distance between right md the next element in Yl. Update right to the next element in Yr. Compute 
distance between left and right. 

The BASIC-2.S4 algorithm is correct and, therefore, the Basic-2 algorithm is also correct, for any 
Minkowski distance with I. 



4 Empirical Time Analysis 

The Basic-2 algorithm has been implemented and tested with randomly generated inputs, starting 
with 125 thousand points and doubling the input size until 16 million points. For each input 
size, 50 independent executions were performed. For each specific generated input, we applied 



the Basic-2 algorithm as well as the standard divide-and-conquer algorithm, described in fCormen 
et al.| pool ], which computes seven pairwise comparison for each point in the central slab. We 
refer to the standard algorithm as the Basic-7, following the naming convention presented in [Jiang 
and Gillespie] pOOT] !. 

Both algorithms were implemented in the C programming language. The source code is avail 
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able atjhttp : //w3 . ualg . ptZ-^f lobo/closest-pair/[ 

The recursion stopped for a number of points less or equal to 10. The algorithms were tested 
with Minkowski distances rfi, d2, di^.uxs and doo. For each independent run, we measured the 
execution time of the entire algorithm. The results are summarized in Figure [9j 



Running tine ratio of Basic-2 over Basic-7 
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Figure 9: Running time ratios of Basic-2 over Basic-7 for various Minkowski distances. 



As expected, and in accordance with the results obtained by [Jiang and Gillespie| POOT] , our 



simulation shows that Basic-2 is faster than the standard divide-and-conquer algorithm, Basic-7. 
Although Basic-2. S4 introduces a few extra relational comparisons (see pseudocode in Figure 
[3]), these are negligible compared to the savings that occur due to the reduction in the total number 
of distance function calls. 

The Basic-2 algorithm is about 20% faster than the standard divide and conquer algorithm 
for the larger input sizes, when using the Minkowski distances d\ and doo. The speedup is more 
pronounced for the case when using the Minkowski distance d2, with Basic-2 being nearly 36% 
faster for the larger input sizes. The reason for the greater speedup when using Minkowski distance 
d2 is because the computation of the distance function is more expensive in this case, and therefore, 
savings in distance function calls have a more profound effect on the overall execution time of the 
algorithm. This fact is fully confirmed by the results obtained when using the somewhat exotic 
distance (is. 1415. Due to the non-integer value of p, the cost of computing this kind of distance 
function is highly inflated and so it is possible to observe speedups over 100% for the Basic-2 
algorithm. 
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5 Final Remarks 



In this paper we analyzed the Basic-2 algorithm, which is an optimized version of the Bentley and 
Shamos procedure for the planar case where the computation of the distance cIlr requires only 
two pairwise comparisons per point to find the closest pair within the central slab. The Basic-2 
algorithm was first presented by [Jiang and Gillesple|]2007| |. 

In Corollary [2] we proved that the two comparisons per point are a minimum for the divide-and- 
conquer CLOSEST- Pair algorithm, for all Minkowski distances with \, This result is a direct 
consequence of the strength of the sparsity condition, which is induced over the set of points in the 
plane by the knowledge of and dR. 

We consider that the generalization of Corollary [2] to higher dimensions (in particular to the 3D 
space) is of interest, considering not only possible applications but also the theoretical significance 
of such an achievement. However, we note that, even in the 3D case this procedure may have rather 
less efficiency gains because of the lack of a natural order relation among the points lying in the 
central slab, and the consequent increase in the problem's complexity. 



References 

Jon Louis Bentley. Multidimensional divide-and-conquer. Communications of the ACM, 23(4): 
214-229, 1980. ISSN 0001-0782. 

M. Smid. Closest-point problems in computational geometry. In J.-R. Sack and J. Urrutia, editors. 
Handbook of Computational Geometry, pages 877-935. Elsevier Science, Amsterdam, 2000. 

Michael Ian Shamos. Geometric complexity. In STOC '75: Proceedings of seventh annual ACM 
symposium on Theory of Computing, pages 224-233, New York, NY, USA, 1975. ACM. 

Jon Louis Bentley and Michael Ian Shamos. Divide-and-conquer in multidimensional space. In 
STOC '76: Proceedings of the eighth annual ACM symposium on Theory of Computing, pages 
220-230, New York, NY, USA, 1976. ACM. 

Jon Kleinberg and Eva Tardos. Algorithm Design. Addison- Wesley Longman Publishing Co., Inc., 
Boston, MA, USA, 2005. ISBN 0321295358. 

T. H. Cormen, C. E. Leiserson, R. L. Rivest, and C. Stein. Introduction to Algorithms. MIT Press, 
USA, 2nd edition, 2001. 

Qi Ge, Haitao Wang, and Hong Zhu. An improved algorithm for finding the closest pair of points. 
Journal of Computer Science and Technology, 21(1):27-31, 2006. 

Minghui Jiang and Joel Gillespie. Engineering the divide-and-conquer closest pair algorithm. 
Journal of Computer Science and Technology, 22(4):532-540, 2007. ISSN 1000-9000. 



13 



